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MRI IMAGING WITH A PERL FIELD 

FIELD OF INVENTION 

This invention relates to MRI (magnetic resonance imaging), and in particular, to 
MRI in which image data is encoded in a spatially-periodic encoding field. 

RELATED APPLICATINOS 

This application claims the benefit of the March 27, 2003 priority date of U.S. 
Provisional Application No. 60/458,048, the contents of which are herein incorporated by 
reference. 

BACKGROUND 

In a conventional MRI system, one often acquires image data in a two- 
dimensional variable space called £-space. The present state of the art steps through each 
line of A>space by stepping through a sequence of gradient fields. To acquire the entire 
2D data set, one must acquire a set of lines in &-space. Since there is a different sequence 
of gradient fields for each line of £-space, this can be a time-consuming process. 

It is known that by encoding image data in a spatially-periodic encoding field, one 
can accelerate the image acquisition. Such spatially-periodic fields are sinusoidal, and 
hence periodic, along one spatial dimension, e.g. x, and vary linearly in the orthogonal 
spatial dimension, e.g. y. These fields have come to be known as "PERL" fields 
("PERiodic in x and Linear in y"). 

U.S. Patent 5,365,172 to Hrovat, et al., the contents of which are herein 
incorporated by reference, discloses a method and system for using a PERL field to 
acquire image data and for processing that data to reconstruct an image. However, 
although Hrovat teaches methods for reconstructing the encoded image, those methods 
rely primarily on numerical methods for inverting an integral transform. 

SUMMARY 

The invention is based, in part, on the recognition of an analytic method for 
inverting a PERL transform. This method enables rapid image reconstruction of MRI 
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images, as well as the solution of problems arising in many other applications that utilize 
the PERL transform. 

Systems and methods according to the invention avoid encoding data using the 
conventional &-space in which k x and k y correspond to two spatial dimensions. Instead, 
the systems and methods encode data in a two-dimensional space characterized by k x and 
m, where k x corresponds to one spatial dimension, as it does in traditional MRI , and m 
represents the index, or order, of a Bessel function. 

The PERL signal includes a superposition of spin echoes weighted by Bessel 
functions of varying order, or index. The process of weighting spin echoes with Bessel 
functions of different order is referred to herein as "m-weighting." These w-weighted 
echoes are separated from each other by varying the phase of a spatially periodic field. 
This spatial phase encoding provides information sufficient to separate out the echoes 
superimposed on each other. Each individual spin echo is naturally encoded in k x space. 
Hence, x spatial information can be obtained in the standard way (by Fourier transform) 
after the superposed echoes are separated. To resolved-dependent information however 
requires a transform along the m-direction. This transform is referred to herein as "the 
PERL transform." 

In one aspect, an image-reconstruction method includes receiving a superposition 
of phase- weighted spin echoes. The superposition is representative of a two-dimensional 
spin density that encodes the image. A partial spin density is then recovered from the 
superposition of spin echoes. This partial spin density is a PERL transform of the two- 
dimensional spin density. The two-dimensional spin density is then recovered from the 
partial spin density by analytically evaluating an inverse PERL transform of the partial 
spin density. Finally, the image is reconstructed on the basis of the two-dimensional spin 
density. 

Particular practices of this invention include those in which recovering the two- 
dimensional spin density includes evaluating an integral of the partial spin density 
weighted by a bandwidth-limited eigenfiinction orthogonal to the Bessel functions. One 
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example of such a bandwidth-limited eigenfunction includes an integral of a Chebycheff 
polynomial. A particular example of such an integral is one of the form: 



where co is a spatial angular frequency and k is a wave number. 

The new methods also include of calibrating a PERL coil by generating a Fourier 
encoded image using a traditional phase encode gradient field in the direction of the 
periodic direction of the PERL coil The PERL field is then substituted for the traditional 
readout gradient field. The reconstructed image has features indicative of any 
misalignment between the PERL field and the phase encode gradient field. The feature's 
locations are then adjusted to reduce the misalignment. This can be achieved by, for 
example, physically moving the PERL coil or by altering the amplitudes of currents 
feeding the PERL coil. 

Other aspects of the invention include methods of calibrating a PERL coil by 
generating a Fourier encoded image using a traditional phase encode gradient field that 
varies along a linear direction of the PERL coil. The PERL field is again substituted for 
the traditional readout gradient field. The Fourier reconstructed image has features 
indicative of a mismatch between a magnitude of the PERL field and a magnitude of the 
gradient field. A current flowing on the PERL coil is then adjusted to reduce the 
mismatch. In this way, the phase encode field, which is already calibrated, can be 
compared with the PERL field amplitude and used to calibrate it. 

Another aspect of the invention includes an apparatus for generating a PERL field. 
The apparatus includes a printed-circuit board having first and second layers. A first trace 
pair has first and second conductive traces that are electrically connected. The first and 
second conductive traces define respective first and second waveforms on the first layer, 
with the second waveform being the mirror image of the first waveform. Both the first 
and second waveforms are spatially periodic in a periodic direction. A second trace pair 
is electrically connected with the first trace pair. The second trace pair has third and 
fourth conductive traces, also electrically connected to each other. The third and fourth 
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conductive traces define respective third and fourth waveforms on the second layer. The 
third and fourth waveforms are identical in shape with the first and second waveforms 
respectively. However, the third conductive trace extends further than the first conductive 
trace along a linear direction orthogonal to the periodic direction. Similarly, the fourth 
conductive trace extends further than the third conducting trace along the linear direction. 

In another aspect, the invention includes a substrate having a first surface and a 
second surface opposed to the first surface. A first conductive trace defines a first 
waveform on the first surface and a second conductive trace electrically connected to the 
first conductive trace defines a second waveform on the second surface. Both the first and 
second waveforms are spatially periodic in a periodic direction. The second waveform is 
identical to the first waveform but spatially offset therefrom in the periodic direction by a 
half of a spatial wavelength of the first waveform. 

Unless otherwise defined, all technical and scientific terms used herein have the 
same meaning as commonly understood by one of ordinary skill in the art to which this 
invention belongs. Although methods and materials similar or equivalent to those 
described herein can be used in the practice or testing of the present invention, suitable 
methods and materials are described below. All publications, patent applications, patents, 
and other references mentioned herein are incorporated by reference in their entirety. In 
case of conflict, the present specification, including definitions, will control. In addition, 
the materials, methods, and examples are illustrative only and not intended to be limiting. 

Other features and advantages of the invention will be apparent from the 
following detailed description, and from the claims. 

BRIEF DESCRIPTION OF THE FIGURES 

FIG. 1 is a graphical representation that shows a PERL field; 

FIG. 2 is a graphical representation that shows an example of an alternating 
current geometry for creating a periodic field. 

FIG. 3 is a graphical representation that shows the dependence of the sine function 
on the number of terms of an expansion. 
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FIG. 4 is a graphical representation that shows contours of the surface shown in 

FIG. 3. 

FIG. 5 is a plot comparing the modified Bessel function with the hyperbolic sine 
function. 

FIG. 6 is schematic diagram of a PERL coil. 
FIG. 7 is a diagram of two conductive traces. 

FIG. 8 is a diagram showing the superposition of the traces shown in FIG. 7. 

FIG. 9 is a diagram of the pulse sequence for PERL imaging. 

FIGS. 10 and 1 1 are images used during calibration of the MRI system. 

FIG. 12 is a graphical representation that shows an exemplary superposition of 
spin echoes detected by the MRI system. 

FIG. 13 is a flow-chart of an exemplary image reconstruction method. 

DETAILED DESCRIPTION 

An MRI system in accordance with the invention uses a spatial encoding field of 
the form: 

B p (x 9 y) = e B G p ysm(qx + 0). [1] 

Such a field, which is shown in FIG. 1, is referred to as a "PERL" field because it has a 
periodic variation in one direction (the x direction in this case) and a linear variation in 
the other direction (they direction in this case). These two directions will be referred to 
as the "periodic direction" and the "linear direction." Although the periodic direction in 
this example is in the x direction, it is understood that the method described herein does 
not depend on this being the case. 

In equation [l],G p is the amplitude of the PERL field and X - 2nl q is the 
wavelength of the spatial variation of the field. The direction of the PERL field B p is 
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parallel to B 0 , the direction of which is given by the unit vector e Bo , which may or may 
not be in the periodic direction without affecting the imaging method disclosed herein. 

The periodic component of a PERL field as shown in FIG. 1 can be created by a 
coil having adjacent currents spaced one-half wavelength apart and directed in opposite 
directions. Such a coil is shown schematically in FIG. 2. 

The PERL imaging sequence uses a spin or gradient echo structure in which the 
PERL field is applied as a pre-encode field and a traditional G x linear gradient is used to 
read out the signal. There is no traditional phase encoding gradient. Since slice-selection 
is along z, the signal S(i) is given by: 

S(t) = ^p(x,y)exp{-iy[G x xt-G p yTsm(qx + 0)]}dxdy 9 [2] 

where Tis the time duration of the PERL pre-encode field and the variable t is the readout 
time. The Jacobi-Anger expansion: 

e acos* = f i rjjA)e imB 9 ' [3] 

m=-oo 

permits an interpretation of equation [2] as the decomposition of a plane wave into a 
series of cylindrical waves. The significance of this expansion lies in the separation of A 
and B into separate factors. For PERL imaging, a more useful form of the Jacobi-Anger 
expansion is given by 

e iAsHB + e) = J m (A)e imB , [4] 

m=-oo 

which explicitly represents the contribution of the spatial phase offset, 0, as a simple 
factor. Using the Jacobi-Anger expansion of equation [4], equation [2] becomes 

S(t) = \\p{x>y)J m W p Ty)dy exp{-iyG x x(t ~ mqljG x )}dx. [5] 

m 

which can be written as a weighted sum of spin echoes 
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with each spin echo being given by the Fourier transform of a partial spin-density 
function: 



Jp«WexpHyG^(<-r m )}£fc, [7a] 



and the one dimensional spin-density function being given by integrating a two- 
dimensional spin density function weighted by a Bessel function, an integral that will be 
referred to as a "PERL transform." 



p m (x)=\p{x 9 y)J m {k y y)dy 9 [7b] 



where, 



k y =j<j p T,md [7c] 

K=rG x {t-t m ) [7d] 

The received signal can thus be viewed as a sum of spin echoes s m (t-t m ), each 
multiplied by a phase factor e md . The m /;i spin echo refocuses at a time 
t m - TE(m) = mqly G x9 at which time its phase is zero. Sequential echoes ( Am = 1 ) are 
separated in time by: 

ATE(Am = l) = ^ = ^- [8] 

To summarize, equations [7a-d] state that the PERL pulse sequence creates a time- 
dependent signal consisting of a time series of spontaneously occurring spin echoes. A 
particular spin echo m refocuses at time t m and can be calculated by first integrating the 
spin density with J m (k y y) as a weighting function and then Fourier transforming the result 
with respect to jc. 
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To achieve a desired resolution, the spin echoes must be sampled at times 
separated by a time interval t given by the Nyquist criterion: 

r= ^ L ~- [9] 

If t is greater than the interval between spin echoes, then neighboring echoes will 
overlap and interfere with each other. Thus the ratio xl ME gives the number of 
overlapping echoes N: 

N = T - = _2*_ = J_ rio] 

A7E(Am = l) g x x res x n 



v res 



Note that l/ATE is the frequency bandwidth over one wavelength of the PERL 
field when the readout gradient is applied, whereas 1/ r is the bandwidth over x res . In 
practice, X > x res . Hence, as a practical matter, it is typically necessary to obtain 
additional information to resolve the overlapping echoes. 

The MRI system is operated with overlapping echoes because, as shown in FIG. 2, 
the periodic component of the PERL coil is produced by parallel current wire segments 
spaced a half-wavelength apart along the periodic direction with current in adjacent wires 
being in opposite directions. Although such a coil produces a periodic field in x 9 
Laplace's equation shows that along z, the PERL field B p decays as exp(-z//l) . Thus, to 

achieve appreciable penetration of the PERL field into a patient, a value of X > x res is 

used. 

Resolution of N overlapping echoes is achieved by acquiring N independent 
signals, each with a different value of the spatial phase shift 0, i.e. 
0 n = 2m IN, n = 0,..., jV - 1 . The signal, i.e., equation [6], acquired from each value of 

the spatial phase encoding (SPE) is denoted by S {n) (t) . These signals are then Fourier 
transformed with respect to 0 n to obtain N time dependent signals with echoes spaced 
apart by NATE = t: 
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SJt) = ^S^(t)e i ^=J j s m+Nk (t-tJ.[ll] 

An MRI system that uses PERL imaging collects data more efficiently than one 
that uses conventional Fourier MRI imaging. In a conventional MRI system, the received 
signal can be represented as 

S(k(t))= llp(x 9 y)cxp{ik(t)r)]}dxdy, [12] 

which clearly shows that if k (t) is a linear function of time, only a portion of &-space is 
encoded. In a conventional MRI system, additional £-space trajectories are required. 
These additional trajectories require time-consuming gradient switching or changes in 
field amplitude. Although faster switching can be achieved by more expensive and 
sophisticated hardware, the resulting rapid changes in the field can induce undesirable 
physiological responses in a patient. 

In contrast, a PERL MRI system provides a received signal represented by • 
S(k x (t))= jlp(x,y)exp{i[k y ysm(qx + 0)- k x (t)x]}dxdy. [13] 

Since dk x (t)/dt is constant, data acquisition no longer requires gradient switching, 
as is the case with the conventional MRI system. The only gradient switching in PERL 
MRI occurs between spatial phase encodes. 

The decay of each m-echo in PERL (see equations [7a-d]) depends only on k x 
(namely, readout with G x ). Since the PERL field is only used for pre-encoding, the 
periodic dependence of the PERL field along x does not create any "uniqueness" 
problems. The two-dimensional pre-encode PERL field results in a two dimensional 
encoding of spatial information. The PERL field is a product of the form f{x)g(y) y a form 
that cannot readily be produced from any combination of linear gradients. As seen in 
equation [5] the PERL field amplitude modulates the spin density (through k y in J n {k y y)) 
and creates an ^-dependent phase shift {m qx) for each m-echo. Applying the readout 
gradient along x thus resolves the various w-echoes by their respective time shifts. 
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The PERL signal has an important symmetry property that reduces the number of 
required spatial encodes. This symmetry property is analogous to the symmetry 
relationship that exists for Fourier MRI, S(k x -k y ) = S*{-k x ,k y ). 

Consider the PERL signal as given by equation [2]. If, in addition to the PERL 
encode, a readout pre-encoding gradient (G xp re) is applied for a duration T x , then 

= {jy*(*>j0exp{^ [14] 

where t x = G xpre T x /G x . Under these conditions the quantity (t-t x ) may be negative. 
Negative times correspond to the collection of echoes with negative m indices (see 
equation [6]). Therefore for a particular value of 0, 

S 0 (-O = {J/?(x,7)exp{z^ 

S e *H)= JJ/totf«pM^ [15] 

s,* (-0 = ^(0 

The signal collected for positive times with a spatial phase encode of 0+ a is the complex 
conjugate of the signal collected for negative times with a spatial phase encode value of 
0. The number of spatial encodes performed in the PERL imaging experiment is equal to 
the number of overlapping echoes to be resolved. By adding a pre-encode readout 
gradient, echoes with negative m indices (corresponding to negative time) can be 
collected. This doubles the acquisition window and halves the number of spatial phase 
encodes, thereby halving the number of gradient switches. 

Information from which an image is reconstructed is contained in the signals 
acquired with different values of the spatial phase. Hence, to reconstruct the image, one 
first performs a Fourier transform with respect to spatial phase to separate the superposed 
echoes. Then, one performs an inverse Fourier transform on each individual echo to 
obtain the partial spin-density function. One can then use the partial spin-density function 
to recover the two-dimensional spin-density function from equation [7b] by performing 
an inverse PERL transform. 
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An analytic method for inverting the PERL transform shown in equation [7b] 
requires functions, M m (k y y) (hereafter referred to as the "M-fiinctions") such that: 

p(y) = I i pJ*)M m (k y y) [16] 

m 

where p(y) = p{y) is the reconstructed spin density. The M-functions have the 
following properties: 

-CO 

t,J m (yW m (y') = Ky-/), [18] 

m=-oo 

where h(y-y') is the system function for the PERL transform and &, m ' is the standard 
Kronecker delta function. The function h(y-y') is essentially the point spread function 
(PSF) for the PERL transform. In terms of the M-functions, the inverse PERL transform is 
the calculated expression for p(y) in terms of p m \ 

P(y)= f jPm M m {k y y) = p*h, [19] 

m=-co 

where /?* h is the convolution of the functions p(y) and h{y). 

The inversion of the PERL transform is related to the expansion of a function in 
terms of Bessel functions of varying index. 

f(kx) = "£a„J n (k 0 x). [20] 

n 

This expansion, which is referred to as the Neumann series, differs significantly 
from the more commonly-used Fourier-Bessel series, in which a function is expanded in 
terms of Bessel functions with varying arguments but constant index. 

Suitable M-functions for inverting the PERL transform are given by the integral of 
a Chebycheff function: 
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M (ky) - - — — fcos m cos" 
kit J 1 



\ k J J 



e iay da). 



[22] 



The M-fUnctions satisfy two important relations: the closure relation: 



k{x-y) 2 n tZ, 



[23] 



and the addition theorem: 



M.{k(x + y))= t,J n - m (kx)M m (ky) 



[24] 



For reference, note that Mo(k y y) is (2/7i)*sinc(^ K y) and that the closure relation 
(equation [23]) is a limiting case of the addition theorem (equation [24]), obtained by 
setting n=0. The M- functions can be used to invert the PERL transform shown in equation 
[7b]. 

The M- functions are also bandwidth-limited eigenfunctions. With this in mind, 
consider the expansion of a function, p(y), which is bandwidth-limited: 



p(y) = ^p m M m (ky) 



[25] 



If p(y) is not bandwidth-limited, a convolution with a sine function, sin(ky)/ky 9 will make 
it so. Because the M- functions are orthogonal to the Bessel functions, the coefficients p m 
can be recovered by evaluating the integral: 



[26] 



Hence, the image is constructed by simple matrix multiplication, as shown by 
equation [25], from the coefficients obtained from the PERL imaging experiment. 
Equation [23] also identifies the sine function as the point-spread function of the system. 
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Resolution and Sampling Constraints 

By defining the maximum spatial frequency that is observable in the PERL 
experiment, the parameter A^in equation [7c] governs the spatial resolution. This is unlike 
standard Fourier imaging, in which the spatial resolution is determined by the field-of- 
view in £-space. 

Like Bessel functions, the functions M n {k y y) are spatially localized. Specifically, 
M n (k y y)*0 in the region \y\ < mlk y . Thus, one can expand p(y) defined in the domain 

-L < y< L using only those M n (k y y) for which 0 < m < k y L . In one embodiment, the 
maximum value of m, M max , used in the series expansion of equation [16] is the smallest 
integer M ^ such that > k y L . 

Because the Bessel functions and the M- functions are both spatially localized, 
when an object in a PERL image moves during image acquisition, any motion artifacts 
generated as a result of that motion corrupt only that portion of the image that includes 
the moving object. This is different from conventional MRI methods that rely solely on 
the Fourier transform. In conventional MRI methods, motion of a portion of the image 
can potentially corrupt the entire image, including those portions that are stationary. A 
PERL based MRI system is thus eminently suitable for real-time or near real-time viewing 
of moving images. 

In another embodiment, M max is determined by examining the dependence of the 
sine function, equation [23], on the number of terms in the expansion. This dependence is 
shown graphically, in FIG. 3, where the function 

Kn,y)=- Y. M •"(W Sty) 

m~-n 

is plotted against n andy. Note that beyond the plateau at which /f=1, the h(n,y) function 
drops off rapidly. 

In FIG. 4, a contour plot shows contours of constant h projected on the ny plane. It 
is apparent from FIG. 4 that h(n f y) is relatively linear. The slope of the 0.99 contour is 
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very close to n. For a field-of-view of 2L with N points, M max > xN/2. This suggests a 
necessary oversampling by id2 along the linear direction. 

PERL Field Geometry 

A variety of PERL fields can be used to encode image information. In one 
embodiment, a PERL field of the form: 



is used. Such a PERL field is created by a single-sided coil (for the exponential decay in 
the slice direction, z) or by a double-sided coil (for the hyperbolic sine and cosine decay 
in the slice direction). An advantage of the foregoing PERL field is that the field-of-view 
in the linear direction (y) is not restricted by the periodicity in the periodic direction (x). 
As a result, it is possible to design coils having a large field-of-view along both the linear 
and periodic directions. 

For the foregoing PERL field geometry, multislice imaging is achieved by 
adjusting the amplitude of the PERL field, G p , for each slice to maintain a constant k y . 

Oblique slice imaging is achieved by imaging the projection of the desired 
oblique slice onto the xy plane. Variation of the amplitude of the PERL field across an 
oblique slice is taken into account by letting the wave number in the linear direction be a 
function of the slice coordinate. For example, suppose a z-slice is tilted along;; (rotation 
about x) in a field whose variation in the slice direction (z) is given by cosh(^z). Then, the 
wave number in the linear direction (see equation [7c]) is given by (k y (z)^ / yG p Tcosh(qz)) 
and the reconstruction along the linear direction;; yields an equation of the following 




[27] 



form: 



m 



[30] 



P(y')= jp(y) T.M m (k >m yy m (k,{x)y)4y 



m 



where k ym is the minimum value of k y in the slice's field-of-view. 
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Application of the addition theorem (equation [24]) gives: 



P(y') = jp(y) M 0 (k ym (y-^-y))dy 



[31] 



For oblique slices it is preferable to generate the PERL field with a single side coil 
to insure monotonicity. A single sided coil also provides greater convenience in 
placement of the coil as a surface gradient coil. 

In another embodiment, a double-sided planar coil can be used to generate a PERL 
field of the form: 



Alternatively, a cylindrical coil can be used to generate a PERL field of the form: 



where I\ is a modified Bessel function. 

Both of the foregoing PERL field geometries share the advantage of depending 
only on two directions, and not on a third direction. Such PERL field geometries are thus 
desirable for three-dimensional, or whole body applications. 

One disadvantage of the foregoing PERL field is apparent from FIG. 5, which 
compares the modified Bessel function with the hyperbolic sine function. As shown in 
the figure, the linear region of both of these functions extends to a maximum qy of 
approximately 1. This results in a field-of-view of approximately X/n. To expand the 
field-of-view, it is therefore desirable to use PERL coil designs that create fields with 
longer spatial wavelengths. Alternatively, the field-of-view can be extended by correcting 
the image for gradient nonlinearities using techniques known in conventional MRL 

The choice of positional variables in all the foregoing PERL fields is arbitrary. 
The various variables can be permuted without consequence. For example, as shown, the 



B z (x, y,z) = G sinh(#y) sin(^x) 



[28] 




[29] 
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periodic direction is the x direction. However, the periodic direction can also be the y 
direction or the z direction. 

In one embodiment, a three-dimensional image is obtained by phase encoding the 
slice direction and performing a two-dimensional PERL experiment at each such phase 
encode. The resulting phase-encoded two-dimensional images are then Fourier 
transformed along the slice direction to produce a three-dimensional image. The total 
number of data acquisitions for this experiment will be N pe *N spe , where N pe is the number 
of phase encodes along the slice direction and N spe is the number of spatial phase 
encodes. 

In another embodiment, a three-dimensional image is obtained by rotating the 
periodic direction of the PERL field 45°. A phase encoding gradient is then applied along 
the slice direction. For example, with the periodic direction rotated to be parallel to the 
line x+y, and with z being the slice direction, the signal for the n th phase encode (with 
gradient G z of duration T z ) is given by 

S n (t) = 2>""' \\\p{x^z)J m {j<} y Ty)dy ex V {-i 7 G x x{t ~™<1 1 4lyG x )-iyGJ 2 z{n-mql42yG z T z )}dxd 

m 

[32] 

As indicated by the above equation, rotating the periodic direction increases the 

spatial wavelength of the PERL field by a factor of 4l , thereby increasing the field-of- 

view. Alternatively the wavelength can be also be decreased by a factor 4l , thereby 
reducing the field-of-view. 

It is also apparent from the above equation that the m-echo maximums will occur 
on the diagonal of the two dimensional n~t space. This increases the distance between m- 
echoes and thereby reduces the number of spatial phase encodes needed to achieve a 
particular resolution. 

FIG. 6 shows a PERL coil formed on a printed circuit board. The printed circuit 
board has a periodic direction and a linear direction, designated as x andy respectively in 



16 



Attorney Docket: 0531 1-025001 



FIG. 6. Each printed-circuit board is a three-layer structure having a total thickness of 
approximately 1.95 mm. 

Each layer of the printed-circuit board has imprinted thereon first and second 
spatially periodic conductive paths. The conductive paths are 8 millimeters wide and 
made of 0.21 millimeter thick copper. The spatial period of each conductive path is one- 
half of the desired spatial wavelength. In the illustrated example, the desired spatial 
wavelength is 5 centimeters. Hence the spatial period is 2.5 centimeters. 

In practice, edge effects arise in the PERL field because the conductive paths 
extend only a finite distance in the periodic direction. To correct for these edge effects, 
the spatial period can be made to vary. The spatial period may either by lengthened or 
shortened to reduce these edge effects. 

The first and second conductive paths on each layer are mirror images of each 
other and are electrically connected such that the same current flows through both paths. 
Preferably, the first and second conductive paths are in series with each other. This mirror 
image geometry allows the linear component of the field to pass through zero in the 
center of the coil. As a result, the PERL field changes sign as one traverses the coil in the 
linear direction. The three pairs of first and second conductive paths, one pair on each 
layer, are electrically connected by vias shown in FIG. 6. 

Although the conductive traces on each layer have the same spatial period, they 
differ in the extent to which they vary in the linear direction. It is this variation in the 
linear direction that causes a linear variation in the field generated by the coil. In the 
illustrated embodiment, these extents differ by multiples of 3.33 centimeters. As a result, 
the current density is highest at the edges of the printed-circuit board and decreases 
towards the center of the printed-circuit board. Because of its dependence on current 
density, the magnetic field increases in amplitude from a minimum at the center of the 
coil to a maximum at the edge of the coil. 

In the embodiment shown in FIG. 6, there are three layers. However, a variation 
of the field in the linear direction can also be affected by only two layers or by more than 
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three layers. In many cases, it is desirable to have a large number of layers, with trace 
pairs on different layers defining paths that are only slightly different from each other. 
Such a configuration generates a magnetic field that more closely approximates a linear 
field in the linear direction. 

The spatial phase of a PERL field generated by a single coil as shown in FIG. 6 
cannot readily be varied. To enable variation in the spatial phase, it is useful to position 
two printed-circuit boards, each configured as shown in FIG. 6, adjacent to each other and 
to offset them in the periodic direction. By appropriately varying the amplitudes of the 
currents in the two printed-circuit boards, it is possible to generate a PERL field having an 
arbitrary spatial phase. 

In one embodiment, the offset in the periodic direction is selected to be 1/4 of the 
desired spatial wavelength. This results in PERL coils that generate PERL fields in spatial 
phase quadrature. 

Another embodiment of a PERL coil is built by superposing trace pairs in which 
each trace has two constituent conductive paths, as shown in FIG. 7. Each path is 
imprinted on a layer of a printed circuit board. 

A path is made up of transverse segments and connecting segments that connect 
adjacent transverse segments. Current on the transverse segments generates a magnetic 
field in the desired direction. The connecting segments are necessary only to establish 
current continuity between adjacent transverse segments. Current on the connecting 
segments generates a magnetic field that is generally orthogonal to the desired direction. 
However, this current also may also generate a transverse component, i.e., a component 
parallel to B 0 , that tends to distort the field generated by the transverse segments. 

To reduce the distorting effect of current on the connecting segments, the two 
paths are superimposed on each other, forming a trace, as shown in FIG. 8, so that the 
currents on the connecting segments generate magnetic fields tending to cancel each 
other. 
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The two paths are connected at vias between the two layers. The alternating 
current elements that produce the periodic component of the field, i.e., the current on the 
transverse segments, are the same in each layer. However, the currents on the connecting 
wire segments joining the transverse are in opposite directions in the two layers. As a 
result, the net magnetic field component arising from current on these connecting wire 
segments is zero. 

The distribution of conducting traces shown in FIG. 7 thus creates a magnetic 
field that is periodic in the periodic direction without introducing anomalous magnetic 
fields arising from current on the connecting wire segments. 

A trace pair made up of traces as shown in FIG. 8 does not, by itself, generate the 
linearly varying component of a PERL field. The linear portion of the PERL field is 
created by including several additional layers of trace pairs made of traces as shown in 
FIG. 8, with each trace having transverse segments of different extent in the linear 
direction. These additional layers are included for each of the two initial layers with 
oppositely directed connecting wire segments. Thus the method for creating the linear 
portion of the field is identical to that described in connection with FIG. 6. 

To achieve a PERL field with arbitrary spatial phase, two coils (referred to as the 
sine and cosine coils) separated by a quarter wavelength can be deployed, as discussed in 
connection with the embodiment of FIG. 6. In some embodiments the sine and cosine 
coils are on different printed circuit boards. However, in other embodiments, the sine and 
cosine coils are on different layers of the same printed circuit board. 

The basic PERL imaging pulse sequence is shown schematically in FIG. 9. Notice 
the existence of a traditional readout pre-encode gradient pulse that forces the m=0 PERL 
echo to occur at the spin echo time, i.e. at the center of the readout gradient. In effect, 
PERL imaging uses two standard linear gradients: one for slice selection and one for 
readout. The conventional phase encode gradient, however, is replaced by the PERL field. 

Calibration of the PERL coil includes matching the characteristics of the field 
produced by the PERL coil with the gradient field produced by a conventional MRI 
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system. This includes physically positioning the PERL coil so that a central zero crossing 
of its periodic field corresponds to the x=0 position of the MRI's G x gradient. It also 
includes selecting a current to pass through the PERL coil so that the amplitude of the 
PERL field thus generated is consistent with the amplitude of the gradient field generated 
by the MRI system. 

To align the position of the PERL field with the gradient field, a standard Fourier 
encoded image is acquired with the conventional phase encode gradient varying along the 
periodic direction, x. However, the readout gradient is not used for recovering the 
encoded image. Instead, the PERL field, which varies linearly in the y direction and 
sinusoidally in the x direction, is used. This results in the formation of an image as shown 
in FIG. 10. The nodes of the waveform in FIG. 10 show the x coordinates of zero 
crossings of the periodic field. To properly position the PERL coil, the x coordinates of 
the nodes in FIG. 10 are measured using the position tool on the scanner, and the PERL 
coil repositioned to locate align a node at jc=0. 

As an alternative to using the position tool and physically moving the PERL coil, 
it is also possible to adjust the x coordinates of the nodes by adjusting the ratio of currents 
in the sine and cosine coils to produce a phase offset. 

Another alternative to using the position tool is to obtain two field maps, one of 
the linear ^-gradient and the other of the PERL field. From the field maps, the zeros of the 
periodic field and the actual zero of the gradient field will be apparent. Again an 
appropriate electronic adjustment to the phase #can be made to zero the offset. 

The foregoing calibration procedure is limited in accuracy to the resolution of a 
display. In some cases, it may be desirable to perform a second calibration step whose 
accuracy is not limited to a display resolution. One such calibration step relies on the 
ability to separate all the spin echoes of a PERL experiment by performing a Fourier 
transform of the spatial phase encoded signals. As noted above, the center of each spin 

echo has the form s m (0)e tme , where s m (0) is a real number, e im0 is a phase factor arising 

from the spatial phase offset, and m is the echo number. The first step is to calculate the 
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phase 0 m at the center of each echo. The 0 m are expected to approximate m0 (for 5^(0) 

positive) or m0 +180° (fors m (0) negative), mapped back into the interval (-180°, 180°). 

Since 0\s small (from the initial calibration), the function m Ovaries slowly with m. This 
enables the 0 m to be fit to the line m0by linear regression. In so doing, it may be 
necessary to add 180° and a multiple of 360° to each 0 m because the corrected phases 0' m 
do not stay in the interval (-180°, 180°), but in fact wind around an arbitrary number of 
times.) 

The error, A# m ,in the measured value 6 m arises in part from noise in the receiver 
coil and from residual overlap between spin echoes. These errors are more severe for 
echoes whose amplitudes are small. Since many more 9 m are available than are needed 
for accurately determining 0 f one approach to reducing error is to ignore the 9 m that 
correspond to echoes whose amplitude is less than a selected threshold. The remaining 
0 m , i.e. those with amplitudes above the threshold, can be used to determine 0. 

The calculated value 0 differs from the actual value by a small error 50. If the 
offset is corrected using the calculated value, the corrected signal 

S(t) = Zs m (t-t m )e imS ° 
m 

will have a residual factor of e im5e . To determine the effect of this residual factor on the 
corrected signal 5(0, it is first necessary to estimate 60. The linear regression gives an 
error 

S0~2yf6A0/(Amf 2 
(assuming a selected threshold equal to the median amplitude, in which case half the 
echoes are used). Thus, m <5# varies by an amount 

Am50~2^f6 A0/(Am) m 
over the full range Am of echoes. If there are Am=m max -m m i n =300-(-300)=600 echoes, 
then for phase measurement errors even as large as A0~1O°, AmS0~ 0.2 A0 y or 2°. 
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To select an appropriate current to pass through the PERL coil, another standard 
Fourier encoded image is acquired, but with the conventional phase encode gradient field 
varying along the linear (y) dimension of the PERL coil. Again, the PERL field is 
sustituted for the conventional read-out gradient. This results in an image as shown in 
FIG. 1 1, the slope of which depends on the ratio of the PERL field amplitude to the phase 
encode gradient amplitude. To determine the current through the PERL coil that produces 
the same bandwidth as that of the phase encode gradient, the current is then adjusted so 
that the angle of the slope is 45 degrees. 

In use, the PERL coil assembly is placed on the patient bed. Because of the 
alternating current structure, there is no net force on the coils, hence no mechanical 
constraints are required to hold it in place. 

In one practice of the invention, a 5-inch receive-only surface coil was placed on 
the bottom inside surface of the PERL coil assembly and a grapefruit was placed inside as 
a phantom. Imaging parameters were: 

• TR = 3000 ms, 

• TE(m=0) = 50 ms, 

• ATE = 96{is, 

• matrix size 96 x 96, 

• slice thickness 2 mm, 

• 16 spatial phase encodes, 

• number of samples per spatial phase encode =2400, 

• maximum number of s m (t -t m ) echoes acquired per excitation (i.e. spatial 
phase encode) m max =192, 

• sampling dwell time =16 microseconds. 
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• readout gradient G x =0.489G/cm, 

• maximum PERL amplitude G p =0. 1 88G/cm 

• PERL encode time duration = 2ms. 

The field-of-view was 300 mm in both the x andy dimensions, resulting in a 
resolution of 3.125 mm x 3.125 mm. The total data acquisition time was the product of 
the number of spatial encodes and TR, which was 48 seconds. 

In the exemplary data shown in FIG. 12, the signal is not refocused at only one 
time point but rather is spread over time. This indicates multiple spin echoes that 
spontaneously refocus at uniformly spaced values of TE. Thus PERL does not require a 
receiver with as large a dynamic range as traditional Fourier methods. 

FIG. 13 shows an overall flow diagram for reconstructing a PERL image in a two- 
dimensional R x R region. 

The first step in the reconstruction method is to generate an R x m max M- function 
table, where m max is the maximum index for M m (k y y). This step occurs only once for a 
given k y value. The M function table is quickly calculated using its Fourier form 
(equation [22]) and FFT (equation [8]). 

The N spatial phase encode signals S M (t) each have a number of sampled points 
N sam = ^max IN > where Nis the number of overlapping echoes. These signals must first 
be phase corrected for any anomalous B 0 field arising from rotational misalignment 
between the PERL coil and the static magnetic field. The anomalous phase is calculated at 
the center of each signal corresponding to the m =0 echo. 

An FFT along the spatial phase encode direction (see equation [11]) increases the 
echo separation spacing to R points, thereby allowing extraction of m max +l echoes (R 
points each) from the initial N x N sam signal array. 
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As seen in equation [6], a first order phase correction with respect to m accounts 
for the spatial phase offset. An FFT of each echo performs the reconstruction along x. 
Finally, the image is reconstructed by matrix multiplication using the M-functions: 

p(x,y) = ^p m (x)M m (k y y). [35] 

m 

Image reconstruction using the analytical form for the M-functions can be 
performed very quickly, thereby allowing timely observation of any data acquired. 

Additional advantages of the method disclosed herein arise from the acquisition 
of data in both the x and y directions in a limited field-of-view without aliasing from 
objects outside that field-of-view. This occurs both because data acquired with the 
readout gradient G x is bandwidth limited and because a signal weighted by a Bessel 
function becomes spatially localized. The methods and systems described herein thus 
enable imaging of limited field-of- views without aliasing. In addition, any motion- 
induced artifacts in field-of-view are spatially localized. They corrupt only a limited 
portion of the image rather than the image as a whole. As a result, the methods and 
systems described herein are well-suited for such applications as, but not limited to, 
coronary artery imaging and all types of image guided therapy. 

OTHER EMBODIMENTS 

It is to be understood that while the invention has been described in conjunction 
with the detailed description thereof, the foregoing description is intended to illustrate 
and not limit the scope of the invention, which is defined by the scope of the appended 
claims. Other aspects, advantages, and modifications are within the scope of the 
following claims. 

The invention has been described in connection with reconstructing an MRI 
image. However, it is to be understood that the invention can be applied towards solving 
problems in many diverse applications where solutions to the form given by equation 
[26] are required. These applications include, without limitation, applications in radar, 
geology, fluid dynamics, acoustics, and optics. 
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Having described the invention, and a preferred embodiment thereof, what we 
as new, and secured by letters patent is: 
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